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Abstract 

We study the effects of friction on the scaling evolution of string networks in 
condensed matter and cosmological contexts. We derive a generalized 'one-scale' 
model with the string correlation length L and velocity v as dynamical variables. 
In non-relativistic systems, we obtain a well-known L oc t 1 / 2 law, showing that 
loop production is important. For electroweak cosmic strings, we show transient 
damped epoch scaling with L oc t ' (or, in the matter era, L oc t 3 / 2 ). A low 
initial density implies an earlier period with L oc t 1 / 2 . For GUT strings, the 
approach to linear scaling L oc t is faster than previously estimated. 



1 Introduction 

The scaling evolution of vortex-string networks has been extensively studied analyt- 
ically in both condensed matter and cosmological settings [1,2,3], but using rather 
different methods. Condensed matter descriptions tend to focus on a coarse-grained 
order parameter </>, providing a low-level picture of defect motion by estimating 
energy dissipation rates. However, high energy physicists take an 'idealized' one- 
dimensional view of string dynamics by integrating out the radial degrees of freedom 
(in the Higgs (j> and other fields) to obtain a low-energy effective action — the Nambu 
action. The resulting equations of motion can then be averaged in order to describe 
the large-scale evolution of a string network in a so-called 'one-scale' model [1]. 
This also involves incorporating energy loss mechanisms, including loop production 
(usually ignored in condensed matter models) and here we add frictional forces due 
to particle-string scattering. Our generalization of this simple 'one-scale' model 
allows us to describe the scaling evolution of a string network in dramatically dif- 
ferent regimes, ranging from a string-forming phase transition in the early universe 
through to a rapid quench in a liquid crystal. (Unless otherwise stated, we use units 
in which Ti = c = Jzb = 1.) 
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2 The evolution equations 



The motion of a cosmic string passing through a background radiation fluid is 
retarded by particle scattering. For a string resulting from the breaking of a local 
gauge symmetry the main contribution comes from Aharonov-Bohm scattering, and 
the corresponding frictional force per unit length is [4] 



where £{ is the 'friction lengthscale', v is the string velocity, /i is the string energy 
per unit length, and T is the background temperature. Here, (3 is a numerical factor 
related to the number of particle species interacting with the string (expected to 
be of order unity for minimal-GUT strings). The phenomenological force law for a 
condensed matter system will take a similar form with the same velocity dependence. 

In an FRW universe with line element ds 2 = a 2 (r) (c/r 2 — dx ), the Nambu 
string equations of motion with the frictional force (1) can be written [4,5] 
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where the 'coordinate energy per unit length' is defined by e 2 = x' 2 /(l — x 2 ), with 
dots denoting derivatives with respect to r, which is identified with the worldsheet 
time, and primes for derivatives with respect to <r, a spacelike coordinate labelling 
points along the string. In deriving (2) and the constraint (3), we require the addi- 
tional gauge condition that the string velocity be orthogonal to the string direction, 
that is, x • x' = 0. This proves to be a particularly useful form because dissipa- 
tion is naturally incorporated in the decay of the coordinate energy density e, while 
preserving the gauge conditions. 

The total string energy is given by E = fJ.a(r) J eda, and we define the average 
rms string velocity v as v 2 = (x 2 ) = J x. 2 eda/ J eda. Then the total string energy 
density p oc E / a 3 will obey the following evolution equation (in terms of physical 
time t, dt/dr = a): 



(4) 



The one-scale model 

Our aim is to study the evolution of the long-string network on the assumption that 
it can be characterized by a single lengthscale L. We therefore draw a distinction 
between long strings and small loops according to whether they are, respectively, 
larger or smaller than L. For Brownian long strings, we can define the 'correlation 
length' L in terms of the network density as p^ = /j/L 2 . 

Following Kibble [1], the probability of a segment of length / moving with a ve- 
locity v intercommuting within a time 8t is approximately IvSt/ L 2 . Consistent with 
our scaling assumption, we then assume that the probability of such an intersection 
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creating a loop of length in the range I to I + dl will be given by a scale-invariant 
function w (l/L). The rate of energy loss into loops is then given by 



dpoo\ v f ( l\ I dl „ poo 

Poo I w — — — = cv— , (5) 



dt J. , V J \LLL L 

/ to loops - J \ / 

where the loop 'chopping' efficiency c is assumed to be constant. (Note that v 
now denotes the average velocity of the long strings only; in previous analyses 
without friction, this velocity was also assumed to be constant and absorbed into 
the definition of c.) Subtracting the loop energy losses we easily obtain the overall 
evolution equation for the characteristic lengthscale L, 

2^-=2HL(l + v 2 ) + ^ + cv. (6) 
at tf 



Velocity evolution 

We now consider the evolution of the average long-string velocity v. A non-relativ- 
istic equation can be easily obtained: it is just Newton's law 

This merely states that curvature accelerates the strings while damping (friction and 
expansion) slows them down. A relativistic generalization of the velocity evolution 
equation (7) can be obtained more rigorously by differentiating the definition of v 2 : 



dt v 



k ( 1 

l- v { 2H + i; 



This involves two assumptions. Firstly, to obtain the damping term we have taken 
(x 4 ) = (x 2 ) 2 ; the difference between the two is a second-order term. Secondly, in 
the curvature term, we have introduced L via the definition of the curvature radius 



vector, 

a(r) „ <i 2 x . . 

where u is a unit vector and s is the physical length along the string. The dimen- 
sionless parameter k is defined by 

((1 -x 2 )(x-u)) =kv(l-v 2 ). (10) 

Clearly k will be affected by the presence of small-scale structure on strings: on a 
perfectly smooth string, u and x will be parallel so k = 1 (up to a second-order 
term as above); however this need not be so for a wiggly string. We shall have more 
to say about k below. 

Equations (6) and (8) form the basis of our generalized 'one-scale' model, which 
we will now proceed to apply in several different contexts. We note that the velocity- 
independent 'one-scale' model (6) has proved to be successful in describing the large- 
scale properties of cosmic string networks in numerical simulations. Any deficiencies 



3 



v vs. t/lf _ L/lf vs. t/lf 




0.5 1 1.5 2 o 0.5 1 1.5 2 



Figure 1 : Log-log plots of the evolution of L and V in condensed matter systems for C = 
(solid lines) and C = 1 (dotted lines), with k = 1 and initial conditions _Lj = ctj = £f, 
Vi = 0.1. The relaxation to scaling is clear. 

seem to be associated with the emergence of significant small-scale structure, that 
is, propagating kinks and wiggles on scales well below L. In friction-dominated 
regimes, therefore, we should anticipate improved quantitative agreement because 
of the suppression of this substructure. 

3 Vortices in condensed matter systems 

As a simple application, we drop the time-dependence of the damping terms in 
(6) and (8), that is, we set H = and £{ = const. This should be appropriate 
for describing vortex-string network evolution in condensed matter systems. The 
evolution equations become (reintroducing the characteristic speed c): 




This yields the following late-time asymptotic behaviour: 

The asymptotic ratio of the loop formation and friction terms is c/k; as expected, in- 
creasing c (or including loop losses in the first place) leads to a lower network scaling 
density and a smaller average velocity v. Fig. 1 illustrates the relaxation to scaling 
for a particular set of initial conditions. The duration and nature of the transient 
regime depends on the initial conditions and the loop 'chopping' efficiency — the 
larger the value of c, the faster the approach to 'scaling'. 

The t 1//2 scaling law for the characteristic lengthscale is a well-known result 
in the theory of phase ordering with a non-conserved order parameter, supported 
by simulations and experiment (see [3] for a recent review). In that context, it is 
obtainde by rather different methods; in particular it should be noted that loop 
formation is not explicitly included. 
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4 Strings in the early universe 

In the early universe the friction lengthscale increases with time (as t 3//2 in the 
radiation era), so friction will only be important at 'early times' [6] . Let T c be 
the temperature of the string-forming phase transition. In a weakly- coupled Higgs 
model G/j = (T c /mpi) 2 , the time of string formation is t c = (fG/j)~ 1 mp^, where 

1 /2 

/ = 47r (7rA/"/45) . Defining t* as the time at which the two damping terms in 
(2,3) have equal magnitude we find 



l /zl\ 3 / 2 R _1 /+ X 1 / 2 

' R*d.), U=l±) TTTT^fMl 11 (Mat.). (13) 



P(G/j) 2 V v ' V37 PI 2 Gii \t Pl 

String dynamics is friction-dominated from t c until £*, after which motion will 
become relativistic or 'free'. A simple heuristic argument due to Kibble (see, for 
example [7]) suggests that in the damped phase the correlation length will scale as 
Loct 5 / 4 . 

As for initial conditions, we assume (following [7]) that in < Li < t{ ~ t c which, 
in terms of our variables, is <r -1 < Li < 1, where a = (3 / / (G/j) 1 ^ 2 and L = L/t c . 
These two extreme limits could respectively correspond to a rapid second-order 
phase transition and a slow first-order transition . Note that using the definitions 
of /3, a and /, one finds that the initial ratio of the string and background densities 
obeys 

^G,<_M <<0.75>, ,14, 



3 V Pb 

(where < u < 1 is a model-dependent parameter). In fact, analysis of the evo- 
lution equations (6,8) shows that for physically reasonable values of c and k these 
bounds hold for all subsequent times. Consequently, cosmic strings can never dom- 
inate the universe. 

We now discuss the different evolution histories of electroweak and GUT strings. 
Electroweak strings 

For strings forming around the electroweak phase transition, we have t c 4 X 
10 32 tpi, t t « 7 x 10 25 t c and a ^ 2 x 10 15 . Note that with = 1 and h ^ 0.7 the 
time of equal matter and radiation densities is t eq ~ 10 22 t c (and the present time 
is to ~ 10 28 t c ), so friction-domination ends well into the matter era. We find that 
the early time evolution of the string network crucially depends on the initial string 
density. The complete evolution from t c through to to is illustrated in fig. 2 for 
several values of Li. 

If the initial string density is low, the correlation length is large, so friction 
dominates over curvature and the strings will move with small velocities. Since both 
the friction and loop formation terms in (62) are velocity-dependent, we find a period 
where strings are conformally stretched, with the dynamical variables behaving as 

--L (-X' 2 v - —- (15) 

tc \t c J crLi t c 

For Li of order unity (close to the upper limit), this stretching regime can last up 
to ten orders of magnitude in time. As the velocity increases, the friction and loop 
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Figure 2: The evolution of electroweak string networks. The plots represent the exponent 
of the power-law dependence of (a) the correlation length L and (b) the average velocity 
V, while (c) shows the actual growth of V and (d) is the ratio of string and background 
densities (the latter two are log-log plots). Solid, dashed and dotted lines correspond to 

different initial conditions with Li = 10 \ Li = 10 4 and Li = 10 8 , respectively; 
the horizontal axis is labelled in orders of magnitude in rescaled time tjt c from the time of 



formation t. 



and also marked are t e q, 



id tr 



formation terms become of the same order of magnitude as the expansion term, and 
there follows a period where we find Kibble's law 



L f2k(c + k) 



1/2 



5/4 



3k 



2(t(c + k) 



1/2 



1/4 



16) 



However, notice that this only holds in the radiation era; in the corresponding 
regime in the matter era, we have L oc t 3//2 and v oc t 1//2 . This latter regime 
is briefly achieved in fig. 2 before friction ceases to dominate the dynamics. The 
network then evolves towards the final scaling regime: 



r _ f9k(k + £) 



1/2 



t 



k 



2(k + c) 



1/2 



(17) 



Electroweak strings only just reach this 'free' scaling regime by the present day. 
Note that if we begin with a sufficiently high string density, then there will be 
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Figure 3: The two 'extreme' cases of approach to scaling for minimal-GUT string net- 
works. Solid and dashed lines correspond to high (Li = 0.05) and low (Li = 0.9) initial 
string density; C and k were chosen to fit radiation-era simulations. Plots represent the 
exponent of the power-law dependence of L and the average velocity V; now t+ ~ 2.93. 

no stretching regime (15) and the L(t) oc t 5//4 regime will begin immediately. A 
shortcoming of our model at extremely high densities is the fact that we do not 
incorporate (otherwise rare) loop reconnections back onto the network. 

GUT strings 

For a GUT model (with Af = 106.75) which produces strings, we have t c sa 3 X 
10 4 tpi, t* sa 855t c , t eq sa 10 50 t c and a sa 29. Fig. 3 illustrates the approach to 
scaling for the two extreme limits for the initial correlation length, 0.035 < Li < 1. 
Since the friction-dominated period between t c and t* is much shorter, the nature 
of the damped evolution becomes much less definite: although we can still identify 
the 'stretching' and 'friction' regimes, there is no longer a clear L oc t 1//2 or I oc t 5//4 
dependence. In particular, the Kibble regime now depends sensitively on the initial 
conditions: the higher the initial string density, the higher the maximum exponent 
in the time-dependence of L. 

The asymptotic 'free' scaling behaviour in the radiation and matter eras is 



e(2 -e) J \e(k + c 

For constant c and k, this model predicts that the ratios of 7 and vo in the matter 
and radiation eras are independent of c and k, being respectively 3/2y/2 and l/y/2 
(whereas simulations, indicate that these ratios are respectively 2 and 0.9 [8]). On 
the other hand, we can reverse the argument and look for the values of c and k that 
match the simulations: 

c r sa0.24, A; r sa0.18, c m sa0.17, k m sa 0.49 . (19) 

If, as we claim, k is related to the small-scale structure (more small-scale structure 
corresponding to a smaller k) these values are qualitatively in agreement with the 
simulations. No doubt we require additional degrees of freedom to treat small-scale 
structure satisfactorily (see, for example, ref. [2]). 
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The main consequence of our analysis is that minimal- GUT strings approach 
the linear scaling regime (and, in particular, become relativistic) faster than pre- 
viously estimated. Another important point concerns the transition between the 
radiation and matter eras which we find is much slower than previously estimated 
(extending for about seven orders of magnitude in time, the deviation from the 
L oc t dependence never being more than 1%). However we must exercise care in 
this 'free' regime given the matching problems at t eq for c and k. 

5 Conclusions 

In conclusion, we have generalized Kibble's one-scale model by including friction 
and treating the average long-string velocity as a dynamical variable. We believe 
this model advances our understanding of the scale-invariant damped evolution 
of string networks, as well as the relaxation to this 'scaling', in both condensed 
matter and cosmological contexts. It should form the basis for further progress 
and, for example, we may be able to gain insight into the effects of string small- 
scale structure via the parameter k (defined in (10)). In a longer more complete 
description, we have also applied this model to the study of global strings or vortices 
(see ref. [9]). 
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